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Abstract 

We study eigenvibrations for inhomogeneous string consisting of two parts with 
strongly contrasting stiffness and mass density. In this work we treat a critical case 
for the high frequency approximations, namely the case when the order of mass 
density inhomogeneity is the same as the order of stiffness inhomogeneity, with 
heavier part being softer. The limit problem for high frequency approximations 
depends nonlinearly on the spectral parameter. The quantization of the spectral 
semiaxies is applied in order to get a close approximations of eigenvalues as well as 
eigenfunctions for the prime problem under perturbation. 
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1 Introduction and problem statement 



Models with high contrasts are widely studied since their unusual properties 
give insight into the behaviour of new meta- and nanomaterials, including 
those which already exist or are reachable nowadays via modern technolo- 
gies. The corresponding mathematical problems often cause computational 
difficulties and require new methods of numerical approximation. A system 
under consideration possessing two components with double high contrasts, 
both in stiffness and mass density, expresses two distinguishing cases of the 
limit eigenvibration behaviour for each of low and high frequency levels. The 
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description of such systems should not be restricted to the construction of 
classical number-by-number eigenfunction asymptotics, which are called low 
frequency approximations. They only ensure close approximations to several 
eigenfunctions corresponding to the bottom of the spectrum. For more pre- 
cise eigenfunction description in the upper part of the spectrum the classical 
approach sets the requirement for e to be negligibly small. Nevertheless, in 
actual physical models the parameter e, denoting the ratio of inhomogeneity 
for a certain physical characteristic, is often small but fixed. Then describing 
actual vibrating systems, a problem of adequate approximation to eigenfunc- 
tions with large numbers arises. In order to solve the problem we propose a 
new asymptotics, being called high frequency approximations, and compare 
them with the classical ones. The high frequency approximations quite pre- 
cisely describe eigenvibrations for which low frequency approximations are not 
precise enough. 

Methods and results. Starting from an operator with a discrete spectrum a 
classical spectral analysis provides the discreteness of low frequency limits for 
eigenelements of the system with high contrasts. Nevertheless, the standard 
approach misses a certain important characteristic, because the completeness 
of eigenfunction system is lost in the limit. Accomplishing the investigation 
and filling up the gaps in the limit behaviour description we construct and 
justify high frequency approximations to the eigenfunctions. The quantization 
conditions play a vital part in the asymptotics providing an e-network on the 
spectral axis in the range of approximation. Therefore even the leading terms 
in the spectral approximations change along with e. Thus we obtain a quite 
precise approximations to eigenelements of the prime problem with a fixed 
small e. Comparing to the previous study of the stiff problems [5], where the 
leading terms of high frequency approximations are independent of e and the 
quantization provides a right choice of the correctors, in the present problem 
the quantization conditions, arising in particular from matching WKB and 
power series expansions, come along with the choice of the leading terms. 

The preliminary results on the limit behaviour of the system under consid- 
eration have been discussed in pQ. The question of asymptotic description of 
low and high frequency eigenvibrations originates in work [2] arising again in 
[3-6] for problems with perturbations of the stiffness only. Elastic problems 
with perturbations of stiffness and mass density, either with other geometries 
or at different perturbation rates, have been studied in [6-12]. 

Problem statement. Let a stiff and relatively light part of the string, which 
occupies an interval (a, 0), be complemented by a flexible and heavy body 
part occupying (0, b) with a < < b. We consider a stiffness coefficient being 
k(x) on (a, 0) and e>c{x) on (0, b), and mass density being er(x) on (a, 0) and 
p(x) on (0,6), with all functions being positive and smooth in [a, 0] and [0,6] 
respectively. We assume that eigenvibrations of the string are described by 
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the self-adjoint eigenvalue problem 



(k(x)u' £ )' + e\ £ r(x)u £ = 0, ie(a,0), u £ (a) = 0, (1) 
e(x(x)u'J + X £ p{x)u e = 0, xg(0,6), u e {b) = 0, (2) 
u e (-0) = u e (+0), fc(0)<(-0) = ex(0)<(+0). (3) 

We investigate the question how the eigenvibrations of the media, namely 
eigenvalues X £ and eigenfunctions it E , change if the parameter e tends to 0. 
More precisely, we look for the good approximations of A £ and u e as e — > 0. 



2 Low frequency approximations 

It is well-known that for each fixed e > the spectrum of problem ([I])-© 
is real and discrete, consisting of simple eigenvalues that form a sequence 
< Af < A| < ■ • • < \ £ n < ■ ■ ■ — > oo as n — > oo. The corresponding eigen- 
functions {u e _ n }^ =1 form a basis in L 2 (a, 6). Moreover, for each number n the 
eigenvalue branch X £ n is a continuous function of e G (0, 1) such that \ £ n < c n e 
with positive constant c n independent of e, which follows from the mini-max 
principle since quadratic forms are continuously depending on e [15] . 

Studying the asymptotic behaviour as e — > of each eigenvalue branch A £ 
with fixed number n and corresponding eigenfunctions u EjTlf , we immediately 
have the convergence £~ 1 A £ — > A„ and u e ^ n — > C/ n , where C/ n = in (a, 0) and 
C/ n in (0, b) coincides with an eigenfunction u\ of the limit problem ([6]) for the 
eigenvalue A n . 

We look for the approximations of eigenvalues and eigenfunctions in the form 

X £ n ~ ejJL n + e 2 v n + . . . , u £ , n ~ u n (x) + ew n (x) + . . . , x G (a, b). (4) 

Constructing standardly the asymptotic expansions we first define the leading 
terms, which satisfy the problem 

for xe(a,0): (k(x)u'J = 0, u n (a) = 0, <(-0) = 0. (5) 

Hence u n = on (a, 0) and therefore 

for are (0,6): (x(x)u' n )' = -fi n p(x)u n , u n (+0) = u n {b) = 0. (6) 

Since we are looking for the eigenfunction approximations, which are sup- 
posed to be different from zero, the limit \i n has to be an eigenvalue with 
corresponding eigenfunction u n of problem (JSJ). 

Let us fix an eigenvalue fi n of (jSJ), and corresponding eigenfunction u n such 



3 



that Jq pu^dx = 1. Then the next terms of (j4j satisfy the problem 

for i6(a,0): (k(x)w'J = 0, w n (a) = 0, (kw' n )(-0) = (x<)(+0). (7) 

Therefore on (a, 0) we have w n = (>cu' n )(+0) J* k~ 1 (t)dt, and 

for x £ (0, 6) : (x(x)w' n )' + p n p(x)w n = -u n p(x)u n , (8) 
w n (b) = 0, w n (+0) = v n (-0). (9) 

The solvability of (jSJ), dSJ) along with ([7]) and the normalization of u n implies 

v n = -(kw n w' n )(-0) = - [°k(w' n ) 2 dx. (10) 

J a 

Then on the interval (0, b) we fix a unique w n such that f b pu n w n dx = 0. 

Justification of low frequency approximations. We use the same letter 
/ both for a function defined on the interval (a, b) and a vector (/_,/ + ), 
where /_, /+ are the restrictions of / to (a, 0) and (0,6) respectively. Let 
£ be the Hilbert space Lf(a, 0) x L^(0, 6) with the scalar product (u,v)c = 

/° ru-V-dx + Jq pu + v + dx and norm = (u,u)]/ 2 , where it = («_,«+). Let 
us introduce the matrix operator ,4 £ in £ 

^ er drr ( ^ drr ) ^ 1 

\ ^ p dx^^ dx^ / 

with the domain 

= {« e C : w_ £ W 2 2 (a,0), u_(a) = 0, u+ e W 2 2 (0, 6), «+(&) = 0, 
u_(0) = «+(0), fc(0)«'_(0) = £x(0)n' + (0)} . 

The *4. £ is a self-adjoint operator with a compact resolvent. The spectrum 
a(A £ ) is the set of all eigenvalues of (jTJ) — (J3J) • 

Let B be a self-adjoint operator in Hilbert space H with a domain T>(B). 
Recall that a pair (p,u) £ R x T>(B) with ||u||# = 1 is a quasimode of the 
operator with an accuracy up to a > if — < a. 

Lemma 1 Suppose that the spectrum of B is discrete. If(p,u) is a quasimode 
of B with accuracy to a, then interval [p — a, p + o~] contains an eigenvalue of 
B. Furthermore, if segment [p — r, p + r], r > 0, contains one and only one 
eigenvalue A of B, then ||w — v\\h < 2t~ 1 <t, where v is an eigenfunction of B 
for the eigenvalue X, \\v\\h = 1. [13] 

Theorem 2 For each n £ N there exists C n > such that the C n e 2 -vicinity 
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of ' e\i n contains exactly one eigenvalue \ £ n of problem 

\K~^n\ <C n e 2 . (11) 

The corresponding normalized eigenfunction u S)Tl satisfies the estimate 
\\u £ . n — u n — ew n \\L 2 ( a ,b) < C n e 2 , with a certain C n > independent of e. 

Proof. We introduce a corrector 0„(x) = a~ 1 w' n (+0)x(x — a) on (a, 0) and 
0n(aO = (0,6) such that = u n + e(w„ + n ) belongs to V(A E ). Let 
A^j = efi n + e 2 u n and = t £ U £ with = HC/^H^ 1 . By the construction 

\\AU E n -KK\\ 2 c< 

K x e\r £ n f\ix n + eu n \ 2 (u' n (+0) 2 + <(+0) 2 ) + K 2 e\r £ n ) 2 u 2 n \\w n \\ 2 Llm , (12) 

with positive constants Kj independent of e and n, and also 

K\ < (1 - e|K + < 1 + C n e (13) 

for e small enough. Therefore a pair and £7^ is a quasimode of «4 e with 
the accuracy up to C n e 2 . By Lemma [H, in C n £: 2 -vicmity of A^ there exists a 
certain eigenvalue A^ of (JTJ) — (J3I) - Additionally, it can be easily shown that the 
eigenvalues converge saving multiplicity, s~ l \ £ n — > \i n . Since the limit problem 
has only simple eigenvalues, in a certain C^e-vicinity of fi n there is no other 
eigenvalues of dH) — (JHJ) that provides ( JTTj) . Applying again Lemma CD finishes 
the proof. □ 

Note that low frequency vibrations vanish in (a, 0) as e — > 0. This naturally 
raises the question on the possibility of constructing other non-trivial on (a, 0) 
approximations of eigenvibrations addressed next. 



3 High frequency approximations 



Considering sufficiently large eigenvalues A^ ~ e 1 (u + eui) 2 with to > 0, we 
look for the asymptotic expansions of eigenfunctions u e>n (x) ~ Y(e,x) with 

v (x) +ev 1 (x) +e 2 v 2 (x), x e (a,0), 

Y(e,x) = { (14) 
(c (x) + eci(x)) siwy £ S(x) + ec 2 (x) cos / y £ S(x), 16 (0,6), 



where Vq is different from zero and 7 e = - + uj\ . The expansion in form (TH 
consists of power series on the interval (a, 0) and two-term short-wave (WKB) 
approximation [15] on (0, 6) since equation (j2J) contains a small parameter 
near the highest derivative. Substituting these expressions into equation and 
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boundary condition ((T]) gives 



(kv' )' + uj 2 rv = 0, v (a) = 0, (15) 
(kv[y + cu 2 rv i = —2uujirv , Vx(a) = 0, (16) 
(fcf^)' + cj 2 rf 2 = —^i^o — 2uuirvi, i>i(a) = 0. (17) 



Next, we substitute F(e, x) into 



e 7 2 (-xS /2 + p) F(e, •) + e 7 ^2x5' C ; + (xS'/co ) cos 7e £ 
+ ecu ( 2xS'c 1 + (xS')'c\ ) cos 7 £ >S 



— e( 2loxS c 2 + c<j(xS' / ) / C2 — (xc )'J cos7 e S' = 0(e ), (18) 

Equating the expressions in the large parentheses to zero we minimize the 
discrepancy in ffTBl . The eikonal equation x5" 2 = p has a solution 



3(x) = t x" 1 / 2 (r)p 1 / 2 (r) dr, x G (0, b). 



X 



Consequently, the transport equation 2xS"c' + (xS")'c = admits a solution 
c(x) = x _1 / 4 (x)p~ 1 / 4 (x) up to a constant multiplier. Therefore c (x) = Pqc(x) 
and Ci(x) = /5ic(x). Introducing h as a unique solution of the problem 

2xS'ti + (xS')'h = (xc')' for x < b, h(b) = 0, 

we set C2(x) = /3ou~ 1 h(x) providing the boundary condition Y(e,b) = is 
satisfied. By construction Y(s, •) formally solves equation ([ID up to the terms 
of order e 3 and equation ([2]) up to the terms of order e 2 . 

We now apply interface conditions ([3]) in order to define parameters u, u 1 , j3 
and Pi. Before that, regularizing the e-dependence of Y(e,+0) we apply the 
restriction 

j+^i) S(0) = 5 + itI, 5 G (-7r/2,7r/2], / G Z. (19) 
Satisfying the interface conditions up to the terms of order e 2 , we set 

«b(0) = (-l) l A)C(0)BiQ«5 

jfc(0K(0) = (-l) / /5 ^ 1 S , (0)x(0)c(0)cos5, 1 J 

ui(0) = (-l) I (/3ic(0)siB < J + yi) 
fc(0K(0) = (-l) z ( / 9icu 1 S / (0)x(0)c(0)cos5 + /3 x(0)(7 2 ), 

where ^ = ^"^(O) cos 5, g 2 = u;iS"(0)c(0) cos 5 + (c'(0) - S"(0)/i(0)) sin 5. 
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Combining (fT5|) and (1201) we obtain that Vq is a solution to the problem 

(kv'Y + uj 2 rv = in (a,0), 

V 7 V ; (22) 

u(a) = 0, k(0)v'(0) sin 5 - wx(0)S"(0)u(0) cos 5 = 0. 

Proposition 3 For even/ c<j > t/iere exists a unique 5 (id) G (— 7r/2,7r/2] 
sitc/i t/iat problem (I22p /ias a nontrivial solution v. 



Proof. If cu 2 is an eigenvalue of the problem (A;?/)' + u 2 rv = 0, f (a) = 0, 
v(0) = 0, we put 5 (id) = 0. Otherwise we consider the eigenvalue problem 

(kv')' + u 2 rv = in (a,0), u(a) = 0, fc(0)v'(0) + //v(0) = (23) 

with respect to the spectral parameter [i. For each u under consideration 
the problem has a unique eigenvalue n(u), which is due to the fact that the 
spectral parameter is missed in equation. Therefore 5 (id) can be found as a 
unique root in (— n/2, 7r/2] of the equation 

u;x(0)|£'(0)|cotano~ = /i(o;). (24) 

Recall that S'(0) < 0. □ 

Fixing an arbitrary id > we also fix w = v(u,x) being corresponding eigen- 
function of non-linear pencil f[2"2"j) with 5 = 5 (id) defined by Proposition [31 Let 
additionally vq be unity normalized in L 2 (a,0). Consequently, (|20|) provides 



We conclude from condition (1221) that the function (3q(uj) is continuous at every 
point u;* for which 5 (a;*) = 0. From (1211) and (j24p we obtain 

(/ct^)' + cu 2 rwi = — 2a;u;irt;o in (a, 0), 
«i(a) = 0, fe(0)ui(0) + /x(w)«i(0) = x(0)/, 

where / = 02sin5(a;) — cjS"(0)oi cos <5(co>) and 5(a>) 7^ 0. The problem admits 
a solution if and only if i> (0)x(0)/ = —2idid\, since /i(tu) is an eigenvalue of 
(I23p . This solvability condition can be derived multiplying the equation by Vq 
and integrating twice by parts. It may be written in the form 

_ (h(0)S\0)-J(0)Bm 2 6(u))v*(u;,0) 
Ul (2u + x(0)S'(0)cos6(u))c(0)sm6(u) 1 ^' * ' 

Thus we get uj\ as a function of id. Additionally, we obtain 
Wl = -A;(0)i7o(0)c 2 (0)(2cj)^ 1 if <?(w) = 0. 
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We now can find vx, which is ambiguously determined. Subordinating it to the 
condition /° rv V\ dx = we fix it uniquely. Then fix is given by (I2~TT) . We fix 
an arbitrary t> 2 being solution of fjl7|) . 

Let us return to condition (1191) . Now it may be considered as the countable 
set of equations for uj: 



Since uj\{uj) can have a vertical asymptote in the interval I = [0, |x(0)|S"(0)|), 
equation (I26p can have roots in I. More subtle analysis shows that for each I 
there always exists a unique root of (125]) in the set [|x(0)|5"(0)|, oo) because 
uj\{uj) — > as uj — > +oo and quantity 5 (a;) is bounded. We consider the roots 
that increase along with /. 

Definition 4 VFe say that uj(1) is an admissible limit frequency for given 
s > and I £ 7* if it is the largest root of (J26]) . 



Let us establish connection between the exact eigenfrequencies yAf and ad- 
dmissible limit frequencies uj(1) in case of constant coefficients. Indeed, for 
k = x = r = g= lwe have ^/A^ = y/eTvn(b — ea) -1 , n G N. Counting 
the admissible frequencies in this case we note that S(x) = b — x and v (x) = 
Cosmuj(x — a) providing, via the proof of Proposition [HI 5(uj) = if uj = nna" 1 
for natural n and cotan 5(uj) = co>~ 1 /i( a ') f° r other cj. Moreover, (1231) yields 
jLt(a;) = — fQ(0)fo(0) _1 = cotan(a;a) gaining 5 (a;) = arctan(tan(a>a)) G 
(— 7r/2,7r/2) or 5 = 7r/2. Observe that c (x) = flo, h(x) = and C2(x) = 
providing cji = 0. Then (|26|) becomes a;6£ _1 = arctan(tan(u;a)) + 7rZ. There- 
fore, ujbe~ x G (— 7r/2 + 7r/, 7r/2 + 7rZ] and tan(u;6£: _1 ) = tan(cja) providing 
uj = enk(l)(b - ea)" 1 for k(l) END Kf with Kf = (z £ (l - 1/2), z £ (l + 1/2)], 
where z e = (1 + ea(b — ea)~ l )~ l . Since z e > 1 and therefore the length 
| iff | is also larger then 1, we have at least one natural k(l) G Kf. Picking 
up the maximal value kY Lax {e) G Kf n N we fix the admissible frequency 



Having exact correspondence for the range of eigenfrequencies and admissible 
frequencies in case of constant coefficients, in general case we further use the 
set of admissible frequencies as the first approximation for the eigenfrequen- 
cies. Let <3> e denote the set of all admissible limit frequencies. The subset $ £ 
of M + is thick enough, the distance between neighboring roots is comparable 
with e. In some sense f[2T)j) could be regarded as a kind of WKB quantization 
condition. The positive spectral ray uj > is covered by the e-net $ £ , for each 
point of which we can construct the asymptotics (ITlj) . For each admissible 
frequency uj G $ e we will denote by Y u (e,x) the corresponding asymptotic 
solution ( JT41 . 




(26) 




uj{1) = enk^^ib - ea)^ 1 . N< 
1<^J. Therefore, ,/aT = ^ 




s 



4 Justification of high frequency approximations 



The function Y u (e,x) can be used to construct a quasimode of the operator 
A £ . Clearly, Y u G £, but Y^ G" V(A e ) because of discontinuity at x = 0. Let 
us introduce functions ( , £1 G C 1 (a, 0) such that Co( a ) = 0, Co(0) = 1 and 
Ci(o) = 0, Ci(0) = 0, Ci(0) = 1- Both functions are extended by zero into 
(0,b). Introducing 

7b(e) = (T w (e, +0) - r w (e, -0)) e" 2 , n (e) = (sY^(e, +0) - 3£(e, -0)) e~\ 

which are bounded in e by construction, we obtain that the function 

Y u (e, ■) = Y u (e, ■) + e 2 (r (£)Co + ri(e)Ci) 

belongs to Setting T w (e, •) = HYJ^ 1 • Y u {e, •) we prove the following 

estimate \\(A £ — e^I)T UJ (e, -)\\ c < Ce 2 by recalling e~ x {oj + sloi) 2 = e^y 2 . 

Proposition 5 The pair (e^ 2 , T u (e, •)) is a quasimode of A £ with accuracy 
to 0(e 2 ) for every admissible frequency ui G <3> e . 

Proposition 6 For the range of numbers n < 9e a ~~ 1 / 2 with arbitrary 9 > 
and < a < 1/2, the eigenvalues satisfy the estimate |A* — sfi n \ < K*e 1+2a . 

Proof. In order to improve (TTTj) we calibrate (I12p . Eigenfrequencis r) n = p}J 2 
and normalized eigenf unctions u n of problem ([6]) cam be represented as [15] 

Vn = j^ + ^ + 0(±) as n-+oo, (27) 
Mn = v /2/ 1 S(0)(xp)- 1 / 4 (sinr /n 5)(l + Ofo" 1 )) on (0,6), (28) 

where ( 128]) is uniform on [0, b] and admits differentiation in x. Then we have 
the approximation of the right-hand side in ([8]), (191 

Un(-0) = WWl + 0(r?~ 1 )), i/ n = -r/ 2 ft 2 ft(l + O^ 1 )), n -> oo, (29) 

with ft = y/2/S(0) (xp)- 4 / 4 (0) and ft = /° k(t)~ l dt. Since there exists the 
fundamental set of solutions corresponding (jSJ) in the form [15] 

y n = (l + 0(^ 1 ))(xp)- 1 / 4 sin^ and g n = (1 + 0(^ 1 ))(xp)~ 1 / 4 costs', 

u> n admits representation «;„ = p n (x)y n (x) + q n (x)g n (x) for certain functions 
p n and g„. Exploring this structure of solution in problem (jSJ), (EJ) we obtain 

w n (a;) = (1 + 0{r]~ 1 ))r} n K{Hp)~ l/A {x) cosr) n S(x), n -> oo, (30) 

with constant if = ftft(xp) 1 / 4 (0) 1 S(0)~ 1 . Then Q3QD and (|2ZD provide 

lk„|| L 2(o,fe) < ifin, \w' n (+ti)\ < K i n - (31) 
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Finally, counting CT-CT in flTJD we obtain 



~ KKWc < K 3 \r s n \e 2 n 2 VT+^ < K^ 2 \ (32) 

for n < 9e a ~ 1/2 and |r*| < 1 - K 5 e a+1/2 , which follows from Then 
application of Lemma [1] finishes the proof. □ 

Theorem 7 Lei 9 > 1, < a < 1/2, < 7 < 1/2 — a. If u — uo(n) is an 
admissible limit frequency from the number range n G [9- 1 e^,9e a - 1 / 2 ] and 
8{u) ^ tt/2 then the eigenvalue A^ and eigenfunction y e ^ n satisfy the estimates 

\\ E n - e~ 1 (uj + auxO)) 2 ! < «i£ 2 , \\y e>n - T w (e, Oil £3(0,6) < a 2 9e 1+1 , 

with positive constants a\, a 2 being independent of e. 

Proof. Let uj £ n = */A^. Proposition Eland (I27p for the given number range yield 

+ ?)+°( £>+T ) and A* = fg + 0(s). (33) 

We now estimate the distance between neighboring eigenvalues of A £ . From 
flUD we have A^ +1 - A* = £7r 2 (2n + 1)^" 2 (0) + 0(e). If n > fl" 1 ^ then 



|A £ n+1 - A £ J > 27r 2 S~ 2 (0)ne + 0(e) > So^V" 7 , (34) 
with constant go being positive and independent of n. By the similar argument, 

- <| > v^(5M - c (^" 2 + ^ +7 )) > v 7 ^ {sk ~ ■ ( 35 ) 

Let for a certain number I the admissible frequency cj* = u(l) minimize the 
difference \\J~Euj e n — (cu* + £a>i(a;*))|, which equals |e7r(n — / — -) + 0(e 1+7 )| by 
(j26j) and (13"3"j) . Note that if |5| < | the latter is minimized only for / = n (for 
sufficiently small e) and then we have 

y/iou E n - K + < \5(u}*)\S~\0)e + 9 2 e 1+ \ (36) 

Suppose that frequency uj^ +1 (the same for u^-J also satisfies the last in- 
equality. Then we obtain the estimate 



i +l -ul\<2\5(u)\S-\Q)V~e + 29 2 e^ 



that contradicts (j3"5|) for e < e because 2|5(o;)| < n. The number e can be 
found from the equation S(Q)e^(29 2 e 1 / 2 + B x ) = tt - 2|<S(u;)|. Hence, cl^ is a 
unique eigenfrequency that satisfies (1561) with uj(n) being the root of (1261) . 

In view of Lemma [1] and Proposition [51 we can improve inequality (1361) to 
|A* — e _1 (o; + eo;i(a;)) 2 | < a\E 2 . Repeated application of Lemma [1] enables 
us to write ||m £) „ — T w (e, OIU < ai2#£ 1+7 , because the spectral gap is of order 
9~ 1 e 1 -~< due to (1331). □ 
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Note that the case S = 7r/2 is not a typical situation. Indeed, in this case the 
admissible frequency u coincides with the eigenfrequency of problem (|22|) with 
Neumann condition v'(0) = 0. In this case we can not establish the vicinity of 
uj that would contain one and only one eigenfrequency J)^. The arguments of 
the last proof show only that there exist not more then two eigenfrequencies 
satisfying f[3"6"j) . If ui e n and uj^ +1 satisfy (J3"Sj) then T w (e, ■) is not yet a good 
approximation to any of u E>n or u £ n+1 . The situation could be improved by 
the next terms of asymptotics but that is out of scope of this paper. 

Numerical example. Let us consider coefficients k = 1 and r = 1 + x 2 on 
the interval (—1,0), x = 1 and p = 1 + x on (0, 1). For the value of small 
parameter e = 0.05, in Fig. 1 we have plotted the eigenfunctions of ([I])-© 
and the leading terms of the low and high frequency approximations given by 
(jl]) and ([UJ) . Let us emphasize that the purpose of high frequency approach is 
in a good approximation of eigenfunctions. In the example only the first few 
eigenfunctions could be well approximated by the low frequency approach. 
Already u e ^ is quite far away from its low frequency limit (see Fig. 1), and 
that can not be improved by the next terms of low frequency asymptotics, 
because the absolute error is large enough. Note that the eigenvalue A§ is still 
far from zero and thus \/A~| can not be treated as a low frequency. In the 
right-hand side plots we observe that the high frequency approximations work 
well for the range of numbers between 5 and 15. We have to mention that the 
proof of Theorem [7] is done by asymptotic methods, so it would be challenge 
to tell in particular examples the exact range of numbers n, for which high 
frequency approximations are valid in case of fixed e. 



We refer to the values of \J\ e n that are calculated with high accuracy as to 
"exact" . The numerical values under discussion are represented in the Table 



n 


5 


10 


15 


exact yA^ 


2.76675 


5.52678 


8.27450 


low freq. approximation ^s/i n 


2.88055 


5.76252 


8.64418 


uj 


0.6270 


1.260 


1.860 




-0.22224 


-0.53779 


-0.02669 


5 


-0.63509 


-1.3217 


-0.03770 


high freq. approximation ^= + y/eui 


2.75433 


5.51464 


8.3122 



As for numerical example we present the low and high frequency approxi- 
mations to eigenfunctions by the leading terms of the expansions only. Thus 
the low frequency approximations yje\i n to eigenfrequencies are given 
in the Table and are accomplished by visualization of eigenfunctions u n of 
problem (jfjj) in the left columns of Fig. 1 (u n is extended by zero to (—1,0)). 
In order to find the admissible frequency ui we need also 5 and u\. In the 
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Low frequency approximations 



High frequency approximations 





Fig. 1. A comparison of the low and high frequency approximations (black plots) 
with the eigenf unctions (grey plots) for u £) §, u £) \q and u £) i5 (from top to bottom). 

vicinity of expected u we create network over uj and for each of this u we find 
S(u) satisfying ( 122]) (up to 1CT 7 ) and then find u\(u) such that ( 125]) has a 
solution. Finally, we find the admissible frequency uo = u(n) giving the best 
approach to (126]) over tabulated u. The high frequency approximations to 
the eigenfunctions we depict from u £ j ~ v(cu,x) for x G (—1,0) and u e j ~ 
Cq(x) sin(- +u>x)S(x) for x G (0,1) with v(u, x) from (122]) . c (x) 



and S(x) = 
ized in L 2 ( 



-1,1). 



x 



sin 5(u) i/l+x 

i 3 ). All depicted eigenfunctions u en are normal- 



Note that the method that is applied for the approximations in 1-dimensional 
case is also applicable in a multidimensional situation. Nevertheless, the jus- 
tification of it requires another technique, which is out of scope of this paper. 
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